Charged pion production in 44RU + 



llRu collisions at 400A and 1528A MeV 
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We present transverse momentum and rapidity spectra of charged pions in central Ru + Ru col- 
lisions at 400^ and 1528^4 MeV. The data exhibit enhanced production at low transverse momenta 
compared to the expectations from the thermal model that includes the decay of A(1232)-resonances 
and thermal pions. Modification of the A-spectral function and the Coulomb interaction are nec- 
essary to describe the detailed shape of the transverse momentum spectra. Within the framework 
of the thermal model, the freeze-out radii of pions are similar at both beam energies. The IQMD 
model reproduces the shapes of the transverse momentum and rapidity spectra of pions, but the 
predicted absolute yields are larger than in the measurements, especially at lower beam energy. 

PACS numbers: 25.75.Dw, 25.75.Ld 



I. INTRODUCTION 

Collision of relativistic heavy ions is a unique method 
to produce a large volume of excited nuclear matter in 
the laboratory at present. At incident beam energies 
near 1 GeV per nucleon nuclear matter can reach its den- 
sity about two to three times higher than normal nuclear 
matter density at the temperature below 100 MeV 
Major motivations for such studies aim to determine the 
nuclear equation of state (EoS) and to study the basic 
properties of quantum chromodynamics (QCD), the the- 
ory of strong interactions. These informations are im- 
portant not only in their own virtue but also for under- 
standing the behavior of astrophysical objects, such as 
neutron stars and supernovae Q]- 

Experiments of this kind started with various ion 
beams at the BEVALAC in 70's Since the early 

90 's the heavy ion synchrotron SIS at GSI-Darmstadt, 
Germany, took over a leading role in relativistic heavy 
ion collisions in the energy range up to 2A GeV. A sce- 
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nario of several complementary experimental setups at 
SIS has allowed to perform thorough investigations of the 
EoS and in- medium properties of hadrons [a B S ^3 ■ 

Several observables, which are accessible in experi- 
ments, have been proposed as sensitive probes to char- 
acterize the properties of hot and dense nuclear matter. 
Am ong the prominent candidates are the collective flow 
[Til [l2j and the particle production . Especially the 
production of pions has been suggested rather early |lj| , 
as their yield can be connected to the temperature of 
the fireball through the nucleon resonances. However, 
later detailed theoretical investigations revealed that pi- 
ons might not be so sensitive to the EoS [TEl IT^ . Nev- 
ertheless the pion production in heavy ion collisions has 
attracted continuous attention because it is the most im- 
portant inelastic channel in nuclear collisions 17, 18] . 
The description of this process is necessary to under- 
stand the whole dynamic evolution of the fireball from 
the early stage to the freeze-out. 

In this paper we present results on the pion produc- 
tion in Ru -I- Ru collisions at 400^ and 1528^ MeV. 
In the analysis we want to test in particular the mod- 
ification of the A-spectral function in heavy ion colli- 
sions ^3 |23| and the Coulomb interaction within 
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the framework of the thermal model. Previously, simi- 
lar ideas have been applied to the pion spectra in Au + 
Au collisions at lA GeV by the KaoS collaboration at 
SIS/GSI [H 113, H?, "22] and at 10.8A GcV by the E877 
collaboration at AGS/BNL 23] . The experimental data 
are also compared to microscopic transport model calcu- 
lations, namely the Isospin Quantum Molecular Dynam- 
ics (IQMD) 24, 25, 26]. Finding the common features 
in the interpretation of experimental data by using these 
two completely different approaches may shed some light 
on the process of particle production in relativistic nu- 
clear collisions. 

In Sec. m the experimental setup and the method for 
the selection of collision centrality are described. The 
main experimental results on the pion spectra as func- 
tions of transverse momentum and rapidity are presented 
in Sec, mil In Secs. lIVI andlVlwe compare the experimen- 
tal data with the thermal model and the IQMD calcula- 
tions, respectively. Finally, conclusions follow in Sec. lVIl 



II. EXPERIMENT 

Collisions of IfRu nuclei with a ^Ru target of 380 
mg/cm^ thickness were studied at 400A and 1528^ MeV 
with the FOPI detector at SIS. The beam intensities 
were typically on the order of 3x10'' ions/sec. A sim- 
ilar amount of events for 'central' and 'medium central' 
conditions was accumulated (some under 'minimum bias' 
conditions). The central trigger required a high multi- 
plicity in the forward Plastic Wall (PLAWA) , which cov- 
ers the laboratory polar angles 9iab = 7° — 27° with full 
azimuthal symmetry, corresponding to about 15 % of the 



total geometric cross section cTgeom- The FOPI detector 
system is described in detail elsewhere |23, H^] • 

For the tracking of charged particles we use the cen- 
tral drift chamber (CDC) which is placed inside a uni- 
form solenoidal magnetic field of 0.6 T. The CDC covers 
diab = 32'^ — 140°. Pions, protons, deuterons and heav- 
ier particles are identified by using the correlation of the 
specific energy loss [dE/dx) and the magnetic rigidity 
(the laboratory momentum p divided by charge) deter- 
mined by the curvature of the tracks. Details of the de- 
tector resolution and performance can be found in Refs. 

[UlllliilsS. 

To illustrate the measured phase space we show the 
CDC acceptance of charged pions at both beam energies 
in Fig. n III this figure, pf^ represents the normalized 
transverse momentum calculated by pt/ {Tn-KPcmlcm), 
where pt is the transverse momentum, m,r is the pion 
mass, and /3cm and are velocity and Lorentz gamma 
factor of the center of mass (cm.), respectively. In addi- 
tion, y'^"^ represents the normalized rapidity yiab/Vcm — 1, 
where yiab and ycm are the pion rapidity in the labora- 
tory frame and the cm. rapidity of the collision system, 
respectively. As a result, the CDC covers more than 90 % 
of the full solid angle when the symmetry of the colliding 
system around midrapidity is utilized. For both beam en- 




-3 -2-10 1-3-2-10 1 2 





1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 

n - 




1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 _ 

-1- : 

n - 




140° ^. 














32° ■ 












■ 




^^^^^ 


" 


3 


-2 -1 


1 


3 


-2 -1 


1 : 

(0) 

y 



to 
oo 

> 



CD 
< 



FIG. 1: Measured raw yields in the plane of normalized trans- 
verse momentum vs. normalized rapidity of vr" (left) and 
(right) in Ru + Ru collisions at 400y4 (top) and 1528A MeV 
(bottom). In all cases the most central 10 % of ageom have 
been selected. Two solid lines in each panel show the geomet- 
rical limits of the CDC. The dashed line in the lower right 
panel represents the upper limit of the pion laboratory mo- 
mentum (800 MeV) in order to separate 7r"'"'s from protons. 
Each successive contour line represents a relative factor of two 
in yields. 



ergies a target absorption eS'ect is visible at small pt val- 
ues near the target rapidity region (—1.2 < j/^"^ < —0.8). 

The collision centrality of each event is determined by 
two methods depending on beam energy. At 400^ MeV 
the variable E^at is used, defined as ratio of total trans- 
verse {E± ) to longitudinal kinetic energy (£J|| ) in the cen- 
ter of mass: 



Er 



(1) 



where i runs over all detected charged particles in the 
CDC and the PLAWA. Previously, it has been demon- 
strated that Erat is a suitable variable for event cen- 
trality, especially in central collisions at beam energies 
< 400^ MeV [sj]. At 1528A MeV the total multiplicity 
seen in the CDC and the PLAWA is used. In this paper 
we select events only for the upper most 380 nib, which 
corresponds to 10 % of the total geometric cross sec- 
tion. Total numbers of analyzed events under these cuts 
are approximately 80000 at 400^ MeV and 240000 at 
1528A MeV. We will adopt the natural units h = c = 1 
in the following. 
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FIG. 2: Invariant spectra of 7r~ (open circles) and vr"*" (solid 
circles) as a function of mt — m,r in Ru + Ru collisions at 
400 A MeV for the most central 10 % of cTgeom- The error bars 
represent the statistical and pt-dependent systematic errors. 
Solid and dashed lines are exponential fit functions {pt > 120 
MeV) to the invariant spectra of tx~ and tt^, respectively. 
The TT^ data and fit functions are scaled by a factor of 0.1 for 
a clearer display. 



III. RESULTS 

Figures 13 and |31 show the experimental charged pion 
spectra in invariant form as a function of the transverse 
mass rat (= + m J ) for several rapidity bins. The 
spectra are corrected for the CDC track reconstruction 
efficiencies which are determined by a GEANT based 
Monte-Carlo simulation 0| taking into account the de- 
tailed response of the FOPI detector; the IQMD model 
was used as event generator |2^. The same method of 
event selection and particle identification was applied in 
the data analysis and in the simulation. The track re- 
construction efficiencies of CDC for tt^ are evaluated by 
taking the ratio of the GEANT output to input spectra. 
This procedure has been applied in the two-dimensional 
space of Pt and y^*^^. On average these reconstruction 
efficiencies are determined to about 62 % and 80 % for 
7r+ and 7r~, respectively, at 400^1 MeV. These values be- 



FIG. 3; Same as Fig. H but for 1528A MeV. The pt- 
dependent systematic errors are negligible at this beam en- 
ergy, so the error bars are dominated by statistics. Solid and 
dashed lines are the sum of two exponential fit functions (pt > 
120 MeV) to the invariant spectra of n" and tt^, respectively. 
The TT^ data and fit functions are scaled by a factor of 0.1 for 
a clearer display. 



come somewhat higher at 1528A MeV, namely ^ 7A % 
for 7r+ and ~ 89 % for tt". Note that the overall track- 
ing efficiency depends on various parameters such as the 
rapidity, pt and the track density. In general, the track- 
ing efficiency for the algorithm used in this analysis is 
worse for smaller momentum. Figures |5] and |31show the 
efficiency corrected spectra. More detailed discussions of 
the systematic errors on the invariant and the rapidity 
spectra will be given later in this section. 

The invariant spectra of pions at 400A MeV exhibit a 
one-slope structure in all rapidity bins (Fig. |2Jl whereas 
clearly two slopes are present at 1528A MeV (Fig. ISJ. 
Consequently, at AQQA MeV, the invariant spectra for 
each rapidity slice are fitted by one exponential function 
as follows: 



1 



2'Kmt dmtdy'^^^ 



C(y(")).exp[ 



(mt - m^) 



(2) 



where C and Tq are the rapidity dependent normalization 
constant and the inverse slope parameter, respectively. 
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TABLE I: Inverse slope parameter To and the measured and the extrapolated numbers of vr^ per event at midrapidity (—0.2 < 
y'"' < 0.0) for the most central 10 % of (Jgeom- The statistical errors are negligible whereas the estimated systematic errors are 
18 and 12 % at 400^ and 1528^4 MeV, respectively. For the later comparison, the pt-integrated results in the same rapidity 
interval from the IQMD calculations are also included. For IQMD, HM and SM represent a hard and a soft EoS, respectively, 
including the momentum dependence of the nucleon interaction. 
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FIG. 4; Rapidity distributions of 7r~ (top) and tt"'" (mid- 
dle) in Ru + Ru collisions at 400A MeV for the most central 
10 % of ageom- Dashed and dotted lines in the two upper 
panels are the calculations for an isotropic thermal source at 
the center of mass and the estimated target /projectile com- 
ponents, respectively. Solid lines represent the sums of two 
contributions. The bottom panel shows the ratio of the rapid- 
ity distributions of 7r~ and tt"*". The error bars in this panel 
consider both statistical and systematic errors, but the dom- 
inant contribution is given by the systematic errors, however. 
The solid horizontal line at 1.32 represents the estimated ratio 
of the isobar model. 



FIG. 5: Same as Fig.|Il but for 1528A MeV. 



At 1528^ MeV, a sum of two exponential functions is 
required to describe spectra. Solid and dashed lines in 
Figs. 121 and |31 represent the corresponding fit functions 
to the n~ and 7r+ invariant spectra, respectively, in the 
range pt > 120 MeV. Note that the one- and two-slope fit 
functions are the simplest phenomenological description 
of the pion data. The fitted Tq values at midrapidity are 
listed in Table U Since the parametrization by Eq. ^ 
describes the pion spectra reasonably well over all mea- 
sured TOt — TO^, they are used to extract the rapidity 
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distributions. In details, we integrate the fitted exponen- 
tial functions bom pt = to oo, hence also accounting for 
the missing pt region in the CDC acceptance (see Fig.^l. 
The extrapolation to the missing low pt region is rather 
significant as shown in Tabled at midrapidity. 

The fully extrapolated rapidity distributions are shown 
in the two top panels of Fig. for 400 A MeV and Fig. El 
for 1528^1 MeV, where the forward cm. spectra are re- 
flected around midrapidity by using the mass symme- 
try of the collision system. The dashed lines in Figs. 01 
and El are the simplest thermal model predictions for an 
isotropic thermal source at the cm. 33]: 



dNi 



th 



oc 



y2 



X cxp 



T coshyc" 
cosh yc ' 



cosh yc 



T 



(3) 



with yc = yiab — Vcm and T = To at midrapidity (higher 
To component at 1528A MeV). Here, T is 35 (44) MeV 
at 400A MeV and 113 (107) MeV at 1528A MeV for tt~ 
(7r+). In Figs. 01 and El the rapidity distributions from 
the isotropic thermal source are normalized at midrapid- 
ity. Obviously, for both beam energies the experimental 
rapidity distributions of tt^ are wider than the isotropic 
thermal source sitting at midrapidity. Also represented 
by dotted lines are the estimated target/projectile com- 
ponents which are merely Gaussian fits to the differences 
between data and the dashed lines. The sums of the two 
contributions are shown by solid lines in Figs. 01 and El 

Our strategy to integrate the exponential fits of the 
invariant spectra for the final dN/dy^^'' distributions has 
been justified by the GEANT based Monte-Carlo simu- 
lation the input rapidity spectra agree nicely with 
the pt-integrated dN/dy'-^^ distributions obtained by the 
exponential fits to the efficiency corrected GEANT out- 
put invariant spectra. Nevertheless the most significant 
source of systematic errors for the dN/dy^'^^ distribution 
is the uncertainty in the estimation of the CDC efficiency 
which relies on the tracking strategy. Such a systematic 
error is significant at 400A MeV, but it becomes much 
smaller at 1528A MeV, however. Furthermore, this un- 
certainty depends on pt. It is negligible in the low pt 
region (< 200 MeV), but increases with transverse mo- 
mentum, reaching about 50 % at pt = 300 MeV in case of 
the AOOA MeV beam energy (the error bars in Figs.|21and 
Elrefiect the systematic as well as statistical errors). As a 
result, the systematic error in the pt-integrated dN/dy'^'^'^ 
distributions which is caused by the uncertainty of the 
CDC track reconstruction is estimated to about 15 % at 
400A MeV and 6 % at 1528A MeV. Other sources of sys- 
tematic errors which are similar at both beam energies 
can be summarized as follows. The systematic error in 
the particle identification (2 %) is determined by chang- 
ing various selection criteria for good event and track 
samples. In addition, different fitting ranges in pt also 
cause maximal 2 % uncertainty in dN/dy^^\ The error 
due to the extrapolation procedure over the complete pt- 
range is estimated to be about 10 %. For this estimation. 



we first use various fit functions, e.g., an exponential fit 
in the invariant and the Boltzmann representations. Fur- 
thermore, we also use the IQMD model; the difference 
between the extrapolated yields by the fit functions to 
the data and the pion yields by IQMD for pt < 120 McV 
are less than 2 % at both beam energies (note that, for 
this test, the IQMD spectra should be normalized to the 
measured pion spectra for pt > 120 MeV because the ab- 
solute yields by IQMD are always larger than the data as 
will be clear in Sec.[VJ. Employing the IQMD model for 
the extrapolation to lower pt can be supported by con- 
firming that the yield ratios of tt~ to 7r+ by IQMD (inde- 
pendent of the EoS) agree with the estimations by the fit 
function to the data within 5 % for pt < 120 MeV at both 
beam energies. Finally, assuming that the sources of var- 
ious systematic errors are incoherent, we calculate overall 
systematic errors for the dN/dy^^'' of both charged pions 
of 18 % and 12 % at 400A and 1528A MeV, respectively, 
by taking a quadratic sum of all contributions. Only the 
statistical errors are shown for the spectra presented in 
this paper, unless explicitly noted differently. 

The pion multiplicities per event at 400A MeV are < 
n^- >= 0.92±0.17(sys) and < >= 0.61 ± O.ll(sys) 



for the most central 10 % of CTo 



The resulting < 



n^- > / < >= 1.51 ± 0.39{sys) agrees with the iso- 
bar model calculation P| , which is shown by the solid line 
in the bottom panel of Fig. 01 within the quoted errors. 
Similarly, the pion multiplicities per event at 1528 A MeV 
are < n^- >= 13.2±1.6(sys), < >= 11.0±1.3(sys), 
yielding < n^r- > / < >— 1.2 ± 0.2{sys), shown in 
the bottom panel of Fig. El again in nice agreement with 
the isobar model result. 



IV. THERMAL MODEL 

In Sec. mil we have seen the indication that the rapid- 
ity distributions of pions in central Ru -I- Ru collisions are 
wider than the isotropic thermal source at midrapidity 
by comparing the data with the simplest thermal model 
formula, Eq. Q '3^ . In order to investigate this observa- 
tion further, a more instructive (but more complicated) 
thermal model is formulated in this section, including the 
decay of A(1232)-resonances to pions explicitely. 

The thermal model adopted in this pap er has been 
described previously in Refs. 0, [23, lU 0. Decay 
pions from A(1232)-resonances, tta, and additional ther- 
mal pions, ttt, are main ingredients in this model. For 
the momentum spectra of A(1232) and ttt at freeze-out, 
we assume an isotropic expanding thermal source at the 
cm. with the freeze-out temperature Tf = 30 31] and 
84 MeV 0] at 400A and 1528A MeV, respectively. Fur- 
thermore, the radial flow velocity /3/ = 0.3 is included 
for both beam energies. We use the following expression, 
which was proposed by Siemens and Rasmussen for the 
first time in late 70 's |35j, for the thermal distributions 
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FIG. 6: Phase shift ^33 (top) and spectral functions of 
the A(1232)-resonance (bottom). In the bottom panel, the 
dashed line represents a normal Breit-Wigner type function 
whereas the solid line represents the thermodynamic spectral 
function including contributions by the tvN interaction pnl ]. 



of A(1232)'s and ttt's: 
1 (PN 



2'Kmt dnitdy'-^^ 



1 lE 

oc E ■ cxp 

V Tf 

Tf . sinh a 

L a 



T 

cosh a], (4) 



where 7 = 1/^/1-/3? and a 



{-/ ■ f3f ■p)/Tf. Here 



E = rrit cosh yc and p = pf + sinh^ yc are the total 
energy and respective momentum of particle in cm. 

Two different mass distributions of the A(1232) reso- 
nance have been considered with and without including 
the ttN interactions in the thermal fireball. Weinhold 
and collaborators have exploited the detailed calculation 
for the thermodynamic potential of a system consisting 
of pions and nucleons 19, 20]. Two A(1232)-spectral 
functions are shown in the bottom panel of Fig. El where 
the dashed and solid lines are for the free and the mod- 
ified spectral functions, respectively. The free spectral 
function An is simply a normal Breit-Wigner shape 



Ar{E,. 



(i;,„-Sfl,)2 + r(s„„)2/4' 



(5) 



where TiEnm) is the energy-dependent width of the 
A(1232) 



r(-Ecm) 



where 



Er 



and 



Er- 





\' Er 




+ 5^\ 








+ 5^) 



fr2 



"N 



^"r 



(6) 



(7) 



(8) 



with TJiff being the nucleon mass. The parameters used 
in this analysis are r° = 120 MeV, Er = 1232 MeV, and 
5 = 300 MeV, following Ref. "2^. 

Now considering the -kN interactions, the spectral 
function is expected to be modified as 



BR{Ecm) — 2 



dE 

cm. 



(9) 



where the phase shift of the A(1232)-resonance (Pas- 
channel), (533, can be deduced by 



tan(533(£;cm) = - 



r(^cm)/2 

Er 



Er 



(10) 



The phase shift factor ^33 is displayed in the upper panel, 
and the resulting A(1232)-spectral function in the bot- 
tom panel of Fig. |H1 It demonstrates a clear difference 
between the two spectral functions Ar and Br, espe- 
cially close to the threshold; the modified function Br is 
shifted to lower masses. Note that Br can be uniquely 
determined by the measured phase shift 533 from Eq. 
in a model-independent way. 

In the present thermal model, the A(1232)-resonances 
in the thermally equilibrated system start to decay at 
the surface of the fireball at freeze-out. The resulting 
invariant spectra of tt^ near cm. (now as a function 
of pt) from the thermal model calculations are shown 
in the left panels of Fig. [7| where the dashed lines are 
obtained by Ar and the solid lines are obtained by Br. 
In the model calculations, thermal weighting factors are 
properly taken into account for both spectral functions. 

In Fig. 13 the ratios of pions from the A(1232)-decay 
to all pions, tta/ (tta + t^t)i are determined by fitting the 
model calculations to the measured pt spectra; ~ 72 % 
at 400A MeV and - 76 % at 1528A MeV. Comparing to 
the previous results, the estimate at 1528A MeV is about 
12 % higher than the number estimated for smaller Ni -|- 
Ni collisions at similar beam energy (1450^ MeV) [33 |. 
This implies that the fraction of projectile and target 
nucleons excited to A(1232)-resonances at freeze-out is 
higher for a larger collision system. But this ratio does 
not show the beam energy dependence between 400^ 
and 1528A MeV. One important aspect is that at 400^ 
MeV the ratio tta /(tta + t^t) is not uniquely determined 
with the transverse momentum spectra alone, as they 
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FIG. 7: (Left) Invariant spectra of tt* at midrapidity as a 
function of pt in comparison with the thermal model calcula- 
tions at 400^ (top) and 1528 A MeV (bottom). The dashed 
line is obtained for the free A(1232)-spectral function whereas 
the solid line is for the modified function (no Coulomb cor- 
rection). The final results of the thermal model for vr" and 
tt"*" are shown by dotted and dash-dotted lines, respectively, 
after the static Coulomb potential effect is taken into account. 
Two thick vertical solid lines represent the low-pt limits for 
the validity of the present thermal model. (Right) Ratios of 
7r~ to tt"*" yields at midrapidity. Two solid lines in each panel 
represent the upper and lower limits utilized in the present 
thermal model calculation. 



can be well described by one exponential function. How- 
ever, by fitting the transverse momentum spectra and 
the yields at midrapidity simultaneously, tt/^/{tt/^ + ttt) 
can be uniquely fixed (see Fig. |S1 and the relevant text 
below) . 

From Fig. [7| it is clear that the free spectral function 
Af( can not reproduce data at low pt region {pt < 250 
MeV), and that the modified spectral function Bji can 
enhance the yields at low pt, which brings the model 
calculation closer to the measurement. The same effect 
has also been observed before at similar beam energies 
between lA and 2 A GeV El IS HI HI . 

The contributions by higher baryonic resonances, e.g., 
A''* (1440), are negligible because of relatively low freeze- 
out temperature. Even in 1.9^ GeV Ni + Ni collisions 
the total contribution by higher mass resonances to total 
pion yield has been estimated to only ~ 5 % hence 
this factor is neglected in the present thermal model anal- 
ysis for the pion yields. Besides the yield estimation. 



higher mass resonances can not explain the differences in 
transverse momentum spectra between the data and the 
model calculations with An in the left panels of Fig. [7| If 
we increase the number of higher mass resonances to com- 
pensate the difference at low pt by N* ^ N + mr chan- 
nels with n > 2, the pion contribution by N* N + tt 
channels also increases accordingly because the branch- 
ing ratios are fixed by the experimental data 38]. But 
the decay momenta of pions in N* N + tt processes 
are usually larger than the ones in the A(1232)-decay so 
that the shape of the spectra at higher pt becomes much 
flatter than the data. Hence we conclude that the modi- 
fication of the A-spectral function is required to describe 
the present pion spectra. 

As next step, we consider the effect of the Coulomb 
potential in order to explain the observed difference be- 
tween the 7r~ and 7r+ spectra in the low-pt region. With 
a static approximation for the Coulomb field '211 US j 
which neglects the time evolution of the fireball, tt^'s at 
freeze-out feel the effective Coulomb potential given by 



Vc 



Rf 



(11) 



where Rf is the radius of the fireball at freeze-out, Z^ff 
the effective charge contained in the fireball. Then, the 
total energy of the emitted pion is modified, depending 
on their charge, in the following way: 



E{p) = E{p,)±Vc. 



(12) 



where pi is the initial pion momentum without Coulomb 
potential. Now the number of particles with momentum 
p, N(p) = d^N/dp^, can be related to N{pi) by the Ja- 
cobian: 



N{p) 



dpi 



dp^ 



Nip^ 



pE{p) 



■N{p,) = C±iV(K) (13) 



by using the identities pdp — EdE and dE{p) — dE{pi), 
where 



C± = Jp^T2E{p)Vc + VS ^(P)T^<^ 



pEip) 



(14) 



for positive (C+) and negative (C~) particles. In or- 
der to determine the strength of the Coulomb poten- 
tial, Vc, the experimental yield ratios of 7r~ to tt"^, i.e., 
N^-(p)/N^+{p), at a given momentum p are fitted by 



y/p^ + 2E{p)Vc + V'S 
- 2E{p)Vc + 
f E{p) + Vc \ N,-{p,) 
\E{p)~Vc) N^^ipi)' 



(15) 



where the normalization constant Cn and the Coulomb 
potential Vc are two free fit parameters to be deter- 
mined by the data. The constant Cp> is responsi- 
ble for the height whereas Vc determines the slope of 
N^- [p) /N^+ (p) at low Pt- The right panels of Fig.[7|show 



8 



the comparisons of the data with thermal model calcu- 
lations at both beam energies, assuming the same mo- 
mentum dependence of tt^ spectra before the Coulomb 
correction, N^-{pi) — N^+{pi) in this case. The val- 
ues of Cr = 0.9 and Vc = 17 ± 6 MeV describe data 
well at 400A MeV, especially in the low-pi region where 
the TT^ spectra differ significantly. Similarly, Cr — 0.9 
and Vc = 16 ± 4 MeV are the best set of parameters at 
1528 A MeV. The agreement between the data and the 
model calculations in A^^-(p)/iV^+(j5) at midrapidity is 
reasonable at both beam energies. Note that the same 
analysis technique as in this paper can reproduce the tt^ 
spectra, measured in Au -I- Au collisions at \A GeV by 
the KaoS collaboration, with Vc = 25 MeV ^f\ . 

It should be emphasized that the dynamical consider- 
ation becomes important only when the velocity of pions 
is smaller than the surface expansion velocity (Jig) of the 
fireball [^^l • Charged pions with the velocity larger than 
/3s see a time-independent Coulomb potential, and the 
present thermal model is applicable only for those pions. 
For the average radial flow velocity /?/ = 0.3, /3s is about 
0.40 (or 0.53) by assuming a linear (or quadratic) flow 
profile as a function of the fireball radius. Then, the crit- 
ical momentum (m^7s/3s) of pion is about 60 (90) MeV 
for the linear (quadratic) flow proflle. These estimated 
limits, shown by thick vertical solid lines in the left pan- 
els of Fig. [7| are close to or slightly lower than the low-pt 
limit of our measured pion spectra, which is about 100 
MeV. 

By using the determined values of Vc, we can generate 
the invariant spectra of tt^ by Eqs. ((T^ and itTH) at 400^ 
and 1528^ MeV. The dotted and dash-dotted lines in the 
left panels of Fig. |3 represent the final invariant spectra 
obtained by the thermal model for tt" and 7r+, respec- 
tively, after the Coulomb correction is included. Good 
agreement between the data and the thermal model cal- 
culations can be found. 

Once having estimated the magnitude of the Coulomb 
potential, one can deduce the source size at thermal 
freeze-out by using Eq. Ullfl . The centrality condition 
used in this analysis corresponds to the average geomet- 
rical impact parameter < bgeom > of ~ 2.3 fm. The 
evaluated number of participant nucleons (Apart) can be 
estimated to about 146 by employing the recipe given in 
Ref. After scaling with Z/A of projectile and target, 
one obtains Z^ff ~ 66.7. The corresponding freeze-out 
radii i?/ = Zefj-e^/Vc are 5.6 ± 2.0 fm at AOOA MeV and 
6.0 ± 1.5 fm at 1528A MeV. These results demonstrate 
that the thermal freeze-out radius is almost independent 
of the beam energy. We also compare the freeze-out radii 
estimated in this paper with two other published values 
that were obtained by using the same method. Having 
applied the present thermal model to the KaoS spectra 
at the SIS, the pion freeze-out radius was estimated as 

6.3 ± 0.5 fm for Au + Au coUisions at lA GeV 
When the same model was applied to the E877 spectra 
at the AGS, the pion freeze-out radius was estimated as 

6.4 ± 4.5 fm for Au -I- Au collisions at 10.8A GeV 



TABLE II: Comparison of the ratio of the estimated pion 
freeze-out radius Rf to the sharp-sphere radius Ro {= 1.2A^^^ 
fm) of the projectile (or target) nucleus at different beam 
energies (Et). 



Et (A GeV) 


System 


Centrality (%) 


Rf/Ro 


0.4 


Ru -1- Ru 


10 


1.0 ± 0.4 


1.0 


Au + Au 


14 


0.9 ± 0.2 


1.5 


Ru + Ru 


10 


1.1 ± 0.3 


10.8 


Au + Au 


4 


0.9 ± 0.6 



Although the error in the AGS result is too big to draw 
any meaningful conclusion, the mean radii do not show 
any beam energy dependence for Au + Au collisions, too. 
It is not possible to compare directly the present results 
with other numbers mentioned above because the system 
size is not the same (the mass number of Au is more than 
twice that of Ru). Therefore, we compare the ratio of the 
estimated freeze-out radius i?/ to the radius of the pro- 
jectile (or target) nucleus. The results are summarized in 
Table Ull the resulting ratios are very similar at different 
beam energies between OAA and 10.8^ GeV despite of a 
small difference in the collision centrality. 

Only the midrapidity pt spectra of pions are tested 
for the present thermal model . At forward and back- 
ward rapidities, the situation is more complicated be- 
cause of the observed incomplete stopping and partial 
transparenc y of pro j ectile and target nuclei at these beam 
energies lill lia li^ . Figure HI shows the measured 
rapidity spectra of tt^'s in comparison with the thermal 
model calculations; the data and the model calculations 
are compared only for pt > 120 MeV due to the uncer- 
tainty of the spectral shape at the low-pt region in the 
present thermal model, which are indicated by thick ver- 
tical solid lines in Fig. [7| This comparison demonstrates 
that the thermal model calculations, assuming isotropic 
midrapidity thermal source, give somewhat narrower ra- 
pidity distributions than the data at both beam energies. 
In order to reproduce the full measured rapidity spectra, 
the longitudinally elongated, rather than isotropic, ther- 
mal source of A(1232)-resonances and pions should be 
assumed. Although this kind of adjustment is, in princi- 
ple, possible by changing the widths of the rapidity dis- 
tributions of A(1232)'s and ttt's, we do not attempt to 
calculate the exact rapidity spectra of pions at present. 



V. COMPARISON WITH IQMD 

The Isospin Quantum Molecular Dynamics (IQMD) 
model is a nonequilibrium transport model which con- 
siders the isospin degree of freedom for the nucleon- 
nucleon cross section and the Coulomb interaction |^|. 
In the framework of the IQMD model, pions are explic- 
itly formed via the decay of the A-resonances, and exper- 
imental cross sections are considered. The IQMD version 
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FIG. 8: Measured rapidity spectra of tt (left) and 
(right) in comparison with the thermal model calculations 
at 400^ (top) and 1528A MeV (bottom). The resuhs from 
the thermal model, shown by the dashed lines, are obtained 
with isotropic thermal source at midrapidity. The modified 
A(1232)-spectral function and Coulomb correction are in- 
cluded. Due to the low-pt limit of applicability for the present 
thermal model (thick vertical solid lines in Fig. |7|l, the data 
and the model calculations are compared only for pt > 120 
MeV. 



utilized in this analysis contains only the lowest three A- 
states. A modified detailed balance treatment, which ac- 
counts for their finite width, is also considered. As a con- 
sequence, the appearing spectral shape of A's is changing 
in the course of the reaction. Therefore, in the charged 
pion spectra, the IQMD model has to show, in princi- 
ple, the effects of the Coulomb interaction in addition to 
other collision effects including the nN interactions. 

IQMD has been successful in reproducing various as- 
pects of heavy ion collisions, including the pion produc- 
tion in Au + Au collisions at lA GeV 26| . In this section, 
we compare the present experimental pion spectra with 
the IQMD calculations in Ru + Ru collisions at 400^1 
and 1528 A MeV. For these calculations, a hard (com- 
pressibility K = 380 MeV) and a soft {K = 200 MeV) 
EoS versions are used, including the momentum depen- 
dence of the nucleon interaction (MDI); they are denoted 
by IQMD(HM) and respective IQMD(SM) in the follow- 
ing. Note that the IQMD events are analyzed with the 
same analysis procedure as the data. As in the data anal- 
ysis, we also use the Erat distribution at AOQA MeV to 
select the most central 10 % of ageom whereas at 1528A 
MeV the total multiplicity measured in the PLAWA and 
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FIG. 9: Comparison of data (circles) with the IQMD(HM) 
calculations (triangles) for the rapidity distributions of tt" 
(top) and tt"'" (middle) in Ru -I- Ru collisions for the most 
central 10 % of CTgeom. Solid lines in these two panels show 
the IQMD(HM) calculations normalized to the total experi- 
mental yields of tt* . The bottom panel shows the ratio of the 
7r~ and vr^ distributions, both for data and model results. 
Results from IQMD(SM) are almost identical to those from 
IQMD(HM) so that we omit in these plots. 



the CDC acceptance is taken. 

Top and middle panels of Figs. and E| show the 
comparisons of the experimental dN/dy^^'' spectra (cir- 
cles) with the IQMD(HM) calculations (open triangles) 
for 7r~ and 7r+, respectively. The IQMD(SM) calcula- 
tions (inverted solid triangle) are almost identical with 
IQMD(HM) at 400A MeV, but they are about 10 % 
larger in yields for both charges at 1528 A MeV. There- 
fore, we include these results only in Fig.^|at the higher 
beam energy. 

In accordance with earlier investigations |4^ . we find 
that the IQMD model produces too many pions in Ru + 
Ru collisions, especially at 400A MeV. In order to com- 
pare the shape of the dTV/d?/^"^ spectra, we need to scale 
the IQMD results by 0.72 and 0.67 for tt" and 7r+, re- 
spectively, at 400A MeV for both HM and SM (solid lines 
in Fig. inj. These scale factors are determined in such a 
way that the integration of the model dN/dy^'^^ yields 
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FIG. 10: Same as Fig.El but for 1528A MeV. Inversed filled 
triangles represent the IQMD(SM) calculation. Dashed lines 
in top two panels show the IQMD(SM) calculations normal- 
ized in such a way that total yields of tt* are the same as 
data. 



FIG. 11: (Left) Invariant spectra of tt* at midrapidity as 
a function of pt in comparison with the IQMD(HM) calcu- 
lations, shown by solid and dashed lines for tt" and tt^, re- 
spectively. (Right) Ratio of tv~ to n'^ yields at midrapidity. 
Crosses are from the IQMD(HM) calculations. Results from 
the IQMD(SM) are very similar to those from the IQMD(HM) 
except about 10 % higher yield at 1528^1 MeV. Top and bot- 
tom panels are for the data and the model calculations at 
400^ and 1528A MeV, respectively. 



the same value as the integrated data distribution. Sim- 
ilarly, we also need to scale the IQMD results by 0.90 for 
TT" and 0.93 for 7r+ at 1528A MeV with HM (solid lines 
in Fig. I10|) . On the other hand, this scaling factor is re- 
duced to 0.82 for 7r~ and 0.84 for tt"*" with the choice of 
SM (dashed lines in Fig.[Tnil. The scaled dN/dy'-°'^ distri- 
butions from the IQMD calculations can nicely reproduce 
the shape of the measured dN/dy^^'' spectra of pions at 
both beam energies. The bottom panels of Figs. [Sjand 
1101 show that the ratios of the dN/dy'^^^ distributions of 
7r~ and tt"*" agree with the data within errors at both 
beam energies. 

The left panels of Fig. display the comparisons of 
data with the model calculations for the invariant spec- 
tra of TT^ at midrapidity. Similar to the experimental 
data, the IQMD(HM) calculations show different shapes 
for TT~ and tt"*" in the low transverse momentum region, 
which can be attributed to the Coulomb interaction. 
These differences can be presented more clearly in the 
ratio N^- {p)/N^+{p) as a function of pt (right panels of 
Fig. Ill|l . The IQMD(HM) calculation agrees very nicely 
with the data. The IQMD(SM) calculation gives very 



similar results on the transverse momentum spectra of 
pions and their ratios except about 10 % higher produc- 
tion than HM at 1528A MeV. 



VI. CONCLUSIONS 

We have presented charged pion spectra in central Ru 
-I- Ru coUisions at 400 A and 1528 A MeV. In both beam 
energies the uppermost 10 % of the geometric cross sec- 
tion have been selected. The measured transverse mo- 
mentum spectra in their invariant form as a function of 
nit can be described very well by one exponential function 
at 400^1 MeV and by a sum of two exponential functions 
at 1528A MeV. 

The results from thermal model have been compared 
to the transverse momentum spectra of pions at midra- 
pidity. Pions from the decay of the thermal A(1232)- 
resonances in addition to thermal pions are main ingre- 
dients of the present thermal model. The momentum 
distributions of A(1232)-resonances and thermal pions 
at freeze-out are assumed to be an expanding isotropic 
source at the center of mass. Pion spectra at low 
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transverse momentum require the modification of the 
A(1232)-spectral function due to the nN interactions at 
both beam energies. The Coulomb potential is necessary 
to explain the different spectral shapes between the tt" 
and 7r+ spectra in the low transverse momentum region. 
With this we find excellent agreement between the data 
and thc> thermal model calculations. Using the estimated 
strength of the Coulomb potential, which is determined 
by the momentum dependent yield ratio of tt" to 7r+, the 
freeze-out radius of the fireball has been deduced to be 
about 6.0 fm at both 400A and 1528^ MeV. The rapidity 
spectra of pions require the longitudinally elongated (not 
isotropic) thermal source of A's and pions. 

Finally, the experimental data are compared to IQMD 
model calculations. Although the absolute yields of 
charged pions in the model are somewhat larger (about 
30 and 10 % at 400^ and 1528A MeV, respectively), the 
shapes of the rapidity and transverse momentum spec- 
tra agree nicely with the data. The IQMD model dis- 
plays that the yield ratios of 7r~ to tt"*" at low transverse 



momenta deviate from unity due to the Coulomb inter- 
action; the agreement between the data and the IQMD 
model is satisfactory. The pion yield is independent of 
the choice of the equation of state at 400^ MeV, whereas 
it is about 10 % higher with a soft than with a hard EoS 
at 1528A MeV . The comparison of the pion transverse 
momentum and rapidity spectra with two completely dif- 
ferent approaches (thermal model and IQMD) renders us 
the importance of the Coulomb interaction and the col- 
lision effects including the TriV interactions in the pion 
spectra at SIS energies. 
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